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Abstract 

Wc apply the Effective Field Theory of Inflation to study the case where the continuous shift sym- 
metry of the Goldstone boson vr is softly broken to a discrete subgroup. This case includes and gen- 
eralizes recently proposed String Theory inspired models of Inflation based on Axion Monodromy. 
The models we study have the property that the 2-point function oscillates as a function of the 
wavenumber, leading to oscillations in the CMB power spectrum. The non-linear realization of time 
diffeomorphisms induces some self-interactions for the Goldstone boson that lead to a peculiar non- 
Gaussianity whose shape oscillates as a function of the wavenumber. We find that in the regime of 
validity of the effective theory, the oscillatory signal contained in the n— point correlation functions, 
with n > 2, is smaller than the one contained in the 2-point function, implying that the signature of 
oscillations, if ever detected, will be easier to find first in the 2-point function, and only then in the 
higher order correlation functions. Still the signal contained in higher-order correlation functions, 
that we study here in generality, could be detected at a subleading level, providing a very compelling 
consistency check for an approximate discrete shift symmetry being realized during inflation. 



1 Introduction 



In the light of the current experimental effort, it is important to explore all the possible 
signatures of Inflation. The Effective Field Theory (EFT) of Inflation [T] is the ideal setup 
for doing this. In fact, as we will briefly review, by realizing that every inflationary model 
spontaneously breaks time-diff.s, it reduces the theory of the fluctuations for the most general 
model of inflation to the one of the Goldstone boson associated with the breaking of time-diff.s. 
As it often occurs with Goldstone bosons, the resulting Lagrangian is highly constrained by 
the symmetries and it allows for a complete phenomenological analysis. The EFT of Inflation 
has been developed in the following papers [IllSlEllllElElEllHllSlIinillllinillSlElIIS]. Many 
new observational signatures, such as the orthogonal shape of the 3-point function [3] , or the 
possibility of having large four-point function without a detectable 3-point function , or the 
possibility to have Lorentz invariant (or conformal invariant) shapes of non-Gaussianities [H], 
have been realized in this setup. The Effective Field Theory for Multifield Inflation has 
been realized in [9], and dissipative effects induced by multiple degrees of freedom have been 
introduced in [T3] . The use of Supersymmetry in the Effective Field Theory of Inflation has 
been introduced in [HI [H]. 

So far, the study of the EFT Lagrangian, both in the context of single clock and of 
multi-field, has been concentrated to the technically natural case where the interactions of 
the Goldstone are protected by an approximate continuos shift symmetry. The purpose of 
this paper is to generalize the study to the case where the continuos shift symmetry of the 
Goldstone boson is broken down to an approximate discrete shift-symmetry subgroup. This 
study is motivated by recent explicit stringy constructions of inflationary models called Axion 
Monodromy |16], where an axion plays the role of the inflaton and inflation occurs as the 
axion slowly rolls down its potential. This potential can be thought of as a standard slow- 
roll inflationary potential, such as for example rn^cf?, with a small superimposed oscillating 
component of the form cos(0/A). These models have the interesting observational consequence 
that all the correlation functions oscillate as a function of the comoving wavenumber. 

While the aformentioned setup relies on the explicit ultraviolet models being considered, 
here we will point out that such kind of self-interactions as well as generalization thereof, can 
be thought simply in terms of effective field theory and symmetries protecting the operators in 
the Lagrangian: an approximate continuous shift symmetry broken to an approximate discrete 
one. This has the advantage of disentangling what we might assume about the UV theory 
and what is actually connected to observations. It further allows us to explore all possible 
signatures in full generality. In the lucky event that we discover some non-Gaussianity or some 
oscillations in the cosmological data, we will therefore be able to connect observations to the 
essential Lagrangian of Inflation, which is the one of the Goldstone boson vr. From there we 
might be able to move up in energy inferring about the UV completion of this Lagrangian. 

In this paper we study the case where the Goldstone boson has an approximate discrete 
shift symmetry. In agreement with former literature, we find that the 2-point function has 
oscillatory features as we change the comoving wavenumber of the modes. Because of the non- 
linear realization of time-diff.s, these oscillatory features translate directly into interactions 
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and therefore non- vanishing higher n-point functions. We systematically study this scenario, 
including limits where the 3-point function is very small and a relatively higher signal is in 
the 4-point function. However, contrary to former literature, we find that, within the validity 
of the effective theory, i. e. without the theory being strongly coupled or without the inclusion 
of new degrees of freedom, and avoiding tuning in the parameters, the signal in higher order 
n-point functions is always smaller than the one in the 2-point function. This implies that 
observational constraints or even detection for these kind of models must come from the 2- 
point function, with only marginal need to analyze higher order ra-point functions, unless 
some (unexpected) degeneracies are present in the 2-point functions. The latter might indeed 
be detectable only at some subleading level, and such a detection would provide an extremely 
compelling consistency check of an approximate discrete shift symmetry for the Goldstone 
boson 71. 



2 Effective Field Theory of Single-Clock Inflation and 
Approximate Shift Symmetry 

In this section we briefly review the effective action for single-clock inflation and we will de- 
scribe how an approximate discrete shift symmetry for the clock field can be implemented 



within this approach. Readers familiar with the EFT of Inflation can skip directly to sec. 2.1 
The effective action was developed in [H |2] and we refer the reader to those papers for a 
detailed explanation. The construction of the effective theory is based on the following con- 
sideration. In a quasi de Sitter background with only one dynamical degree of freedom, there 
is a privileged spatial slicing given by the physical clock which allows us to smoothly connect 
to a decelerated hot Big Bang evolution. The slicing is usually realized by a time evolving 
scalar (pit), but this does not need necessarily to be the case. To describe perturbations 
around this classical solution one can choose a gauge where the privileged slicing coincides 
with surfaces of constant t, for example 6(f){x, t) = 0. In this 'unitary' gauge there are no 
explicit scalar perturbations but only metric fluctuations. As time diff.s have been fixed, the 
graviton now describes one additional degree of freedom: the scalar perturbation has been 
eaten by the metric. One therefore can build the most general effective action with operators 
that are functions of the metric fluctuations and that are invariant under the linearly-realized 
time-dependent spatial diff.s. As usual with effective field theories, this can be done in a 
low energy expansion in fluctuations of the fields and derivatives. We obtain the following 
Lagrangian [I1I2]: 

^E.H. + S.F. = J d\ ^\^Ml,R + Mliif^oo - Ml,{2>H^ + H) + 

+ ^M,(t) V° + If + ^M3(t) V° + If + 
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where we denote by 6K^,y the variation of the extrinsic curvature of constant time surfaces 
with respect to the unperturbed FRW: 6K^i^ = Kf^^ — a^Hh^,^ with h^^, being the induced 
spatial metric, and where M2,3 and Mi,2,3 represent some time-dependent mass scales. 

Let us comment briefly on ([T|. The first term is the Eistein-Hilbert action. The first three 
terms are the only ones that start linearly in the metric fluctuations. The coefficients are 
such that when combined the linear terms in the fluctuations cancel. The action must start 
quadratic in the fluctuations. The terms in the second line start quadratic in the fluctuations 
and have no derivatives. The terms in third line represent higher derivative terms. Dots 
represent operators that start at higher order in the perturbations or in derivatives. This is 
the most general action for single clock inflation [T]. 

The unitary gauge Lagrangian describes three degrees of freedom: the two graviton he- 
licities and a scalar mode. This mode will become explicit after one performs a broken time 
diffeomorphism (Stiickelberg trick) to reintroduce the Goldstone boson which non-linearly 
realizes this symmetry. In analogy with the equivalence theorem for the longitudinal compo- 
nents of a massive gauge boson [T7j, the Goldstone decouples from the two graviton helicities 
at high energies, and the mixing can be neglected. As we will review and explicitly check 
later, it is possible to verify that in most situations of interest this is indeed the case and one 
can neglect the metric fluctuations. 

As anticipated, we reintroduce the Goldstone boson (tt) by performing a broken time-diff., 
calling the parameter of the transformation — tt, and then declaring vr to be a field that under 
time diff.s of the form t ^ t + ^^{x) transforms as 



TT X 



^ n{x{x)) = 7i{x) -^^{x) . (2) 



In this way diff. invariance is restored at all orders. For example the terms containing g^'^ in 
the Lagrangian give rise to the following terms: 

r ^ ^^^^^'^^ ^ r + 2g^'^d,n + id.f. (3) 

We refer to [1] for details. If we are interested just in effects that are not dominated by 
the mixing with gravity, then we can neglect the metric perturbations and just keep the vr 
fluctuations. In this regime, a term of the form g'^^ in the unitary gauge Lagrangian becomes: 

^00 ^ _i_2vr-7V2 + l(5,7r)2 . (4) 

Furthermore, we can assume that the vr has an approximate continuous shift symmetry, which 
becomes exact in the limit when the space time is exactly de Sitter [1]. This allows us to 
neglect terms in tt without a derivative that are generated by the time dependence of the 
coefficients in ([T]). Implementing the above procedure in the Lagrangian of ([T|, we obtain the 
rather simple result: 



S'TT — I d OC 



(5) 
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where for simplicity we have neglected the terms originating from the extrinsic curvature as 
they are usually only important for inflationary models where the space time is very close to 
de-Sitter space [1]. 

We notice that when M2 is different from zero the speed of sound of the fluctuations is 
different from one. We have the following relation: 



1 - cl MlM 



Ml = —^-^. (6) 



At leading order in derivatives, there are two independent cubic self-interactions, 7i{diTTy and 
TT^, which can induce detectable non-Gaussianities in the primordial density perturbations. A 
small speed of sound (i.e. a large M2) forces large self-interactions of the form 7r(c}j7r)^, while 
the coefficient of the operator tt^ is not fixed because it also depends on M3. Cosmo logical 
data can therefore constrain (or measure) the parameters of the above Lagrangian. This 
approach has been recently applied to the WMAP data in [3], giving constraints on M2 and 
M3, as well as on the higher derivative operators that we have omitted in ([s]). This is exactly 
analogous to what happens for data from particle accelerators when the precision electroweak 
tests of the Standard Model are carried out [181 |T9| . 



2.1 Approximate discrete shift symmetry 

When above we obtained eq. (|5]), we explicitly assumed that there was an approximate con- 
tinuous shift symmetry for vr. This allowed us to neglect all the terms that would result from 
the time dependence of the coefficients present in the action, as, when we reinsert tt from 
unitary gauge, any function of time becomes function of vr, explicitly: 

f{t)^ f{t-n)^f{t)-fn + ... . (7) 

We notice that neglecting these terms has nothing to do, at least in principle, with taking the 
decoupling limit. It is simply a technically natural assumption of imposing this symmetry on 
TT. However, it is conceivable that during inflation this symmetry could be broken. This is so 
because inflation needs to be a phase of quasi de-Sitter space once averaged on time scales of 
order H. This is indeed the necessary requirement for solving the horizon problem, and also to 
produce quasi scale-invariant perturbations, at least in each reasonable bin of wavenumbers. 
Scale invariant perturbations are a consequence of time translation during inflation, and it is 
conceivable that time translation could be broken on short time scales and be approximately 
recovered on long time scales, to give rise to experimentally acceptable quasi scale-invariant 
perturbations. On smaller time scales, the space time can be quite different from de Sitter, 
and still give an acceptable model of inflation that predicts scale invariant perturbations. 

Since time translation of the background is mapped into shift symmetry of vr, this discus- 
sion leads us to explore the possibility of relaxing the continous shift symmetry for vr. It is in 
general very hard to protect the lightness of a scalar field, and shift symmetries represents a 
natural way to do this Breaking the shift symmetry in general means that there are strong 

^Another possibility is Supersymmetry, but in general it is not powerful enough in an expanding universe. 
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corrections to the mass of the fields that make its hghtness fine tuned. However, there is 
at least one known exception to this. This is the case when the field has a violently broken 
continuous shift symmetry, which however leaves a softly broken discrete shift symmetry. Let 
us consider for example a scalar field in Minkowski space with a Lagrangian 



This Lagrangian is typical for axions. This theory is non-renormalizable with unitarity bound 
equal to 



In order for the particle to have mass smaller than the unitary bound, we need to impose 
II < F. The parameter /i softly breaks the continuous shift symmetry of (f) leaving out an 
unbroken discrete shift symmetry 



Notice that any interaction breaking the continuos shift symmetry is proportional to /i^. 
Possibly there could be operators of the form ~ ((90)"/-F^"~^ or other operators induced 
by loops of the potential terms that are compatible with the continuos shift symmetry and 
therefore cannot renormahze operators that break it. This means that ^ ^ F is technically 
natural. The symmetry pattern in this case is such that at energies n <^ E <^ F, the 
continuous shift symmetry is approximate, while the discrete shift symmetry is always exact 
down to energies comparable to the mass of the particle /jl'^/F, which is the lowest energy 
scale at which the effective theory makes sense. 

Let us now softly break the discrete shift symmetry. Soft breaking means that the radiative 
corrections induced by the terms breaking the symmetry are small and leave the original 
symmetry to be a good approximate one. Let us call the scale suppressing the operators 
breaking the discrete shift symmetry F^. Examples of such terms could be a term as (p^/F^. 
In this case we need to ensure that the loop-induced potential operators from this term 
are smaller than the mass induced by the oscillatory potential. This implies the constraint 
Fd > F^/iJ^. Another possibility for softly breaking the discrete shift symmetry is by adding 
a potential of the form /x^ cos(0/Frf), which is stable under radiative corrections. Here the 
condition that the original discrete shift symmetry is a good approximate one translates into 
Frf > F for /id ^ or > Ffi^/ fi^ for fid ^ /U. This is a pattern of symmetry breaking when 
a continuos shift symmetry is softly broken to a discrete one that is in turn is softly broken. 

The purpose of this paper is to apply the above pattern of symmetry breaking to the 
Effective Field Theory of Infiation. It is technically natural for the Goldstone boson tt to 
have an approximate continuous shift symmetry. However, as we just discussed, it is also 
technically natural to have a softly broken discrete shift symmetry. The simplest scenario 
in which this can be realized is obtained by setting to zero the coefficient of all higher order 
operators in unitary gauge (we will include them later) and taking H{t) and H{t) to have 




(8) 



(9) 



(j) ^ + 2ttF . 



(10) 
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small oscillating components: 



S = / (fx^/^M^^ 



- (3H\t + tt) + H{t + vr)) + H{t + tt) (-(1 + Tr)^ + {d,nf) 



(11) 



Here we have neglected metric fluctuations. We will come back later to this point arguing 
that they are irrelevant. The original models considered in the literature [23| 12^ [251 l26t [T6l 
where a small sinusoidal term has been added to the slow-roll potential 



S = d^x a/— ^ 



^(90)2-ysr(0)-/cOs(0/F) 



(12) 



are included into this category. Here V^r is the slow roll potential breaking the discrete shift 
symmetry of the inflaton (j). 

The mentioned pattern of symmetry breaking is realized if we imagine that the Hubble 
scale is the superposition of a slowly time-dependent function and a small rapidly oscillating 
function: 

H{t) = H,,it) + i^osc(t) sin(a;t). (13) 

Here Hsr(t) and Hosdt) are slowly time dependent. Their time dependence is parametrized 
by the slow roll parameter e 

which is the parameter controlling the smallness of the breaking of the discrete shift symmetry. 
In order to have inflation, we require ifosc ^ -f^sr- Indeed e controls the breaking of the 
discrete shift symmetry. The continuous shift symmetry is broken by i^osc to a discrete one 
7r{x,t) — )■ 7r(af, t) +2tt/u. We are interested in the regime where the discrete shift symmetry 
is softly broken, which means that the unitarity bound induced by the operators that are 
compatible with the discrete symmetry is smaller than the one induced by the operators that 
break it. If we substitute back H{t + vr) into the action, and Taylor expand to obtain the 
leading operators, we will easily see that the relative weight of operators compatible with 
the discrete symmetry to the ones not respecting it is the ratio of the time derivative of 
the oscillating and non oscillating parts of H. Requiring therefore that the cutoff from the 
shift-symmetry-respecting operators ~ F ^ (— 2ijMpj)^/^/a; is smaller than the unitarity 
bounds induced by the gravity induced interactions (~ Mpi) and by the slow-roll mediated 
interactions, we obtain the condition 

« - I » e^/^ ■ (15) 

It is worth highlighting the hierarchy among the various components of H as the various 
powers of the time derivatives are considered. As we mentioned before, the non-oscillating 
part of H is dominant. This must be the case also for H as otherwise the sign of H would 
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change and the Goldstone boson tt would become a ghost, leading to a catastrophic instability. 
This implies e > a Hose/ H and suggests us to define 

"^^OSC^^ II / \ 

eosc = , ICoscI < 1 • (16) 

The oscillatory term can potentially dominate only starting at the level of H. 



2.2 Two-point function 

Let us now compute the 2-point function. In the decoupling limit where we neglect metric 
fiuctuations, the Taylor expansion of the action in (11) reads 

3 



Sn 



0-{n+l)^n 



[n-l)V 



Tj(n) ri—l- 
-11^ '71 71 



[n 



2)\ 



(17) 



(n-2)! 



(0,71)'') + ^(2M(") - 9n^'))7r" 



where in the third line we have integrated by parts TT^-^TT and = d'^H. We will justify 
later having taken the decoupling limit. The linear equation of motion reads 

n + {3H + H/H)7T - a^TT = , (18) 

where we neglected a term in e = —H/H^. This also will be justified shortly. Introducing 
the conformal time t = / dt/a, and defining x = — /cr, eq. (18) can be written as 

ax - ^^^a,.7rk + TTk = , (19) 

X 

where 5 = H / {2HH). This equation is very similar to the one found in [201 El] for the variable 
S(j) and can be solved in a very similar way. So we delay the derivation to Sec. 4.2.1 The 
important point is that S oscillates with time. As modes start deep inside the horizon and 
redshift up to Hubble, their proper frequency crosses u before Hubble crossing if a > 1. In 
this case the modes go through a resonance and the state becomes different from the Bunch- 
Davies vacuum while still well inside the horizon. In the a ^ 1 case the solution to the linear 
equation can be found in the saddle point approximation and by expanding at linear order in 
Cose- We will shortly see that eosc -C 1 for observational constraint. The effect on the power 
spectrum of C = —H7!- + 0(7r^) is given by 

r2/ 



1 - 



7r\i/2 



. a 



1/2 



sin (Q;ln(2A;/A;^)) 



a > 1, 



(20) 



where t^, is such that k/a{t^,) = H{t^) and we have assumed that the time of the resonance 
and the time of horizon crossing are not too far, so that variations of ifsr , ■ ■ ■ between the 
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two times are negligible. Here and in the rest of the paper is a comoving wavenumber 
representing the physical phase of the oscillating term. The leading effects from the slow 
roll parametes have been included with the standard treatment of evaluating all quantities at 
horizon crossing. This result agrees with the one in [201 EI] |5 

Notice that the effect of oscillations appears as a modification of the tilt: a sort of oscil- 
lating tilt. Constraints on the tilt imply that (see [20] for a detailed analysis) 

,1/2 < 10-2 



~ e 



eosc < e ~ 10- 



< 1 , for a > 1 



(21) 



which justifies our assumption of expanding at linear level in eosc- In the case a < 1, the 
requirement of soft breaking of the discrete shift symmetry implies a > e^^^, which combined 



with (21) leads to the bound 



< e'/' « 1 



(22) 



We can therefore treat eosc as an expansion parameter much smaller than one for all values 
of a. 

In the case a -C 1, the oscillations of the background is slow with respect to Hubble, and 
the modes do not undergo through a resonance before freezing out. In this case the effect of 
the oscillations shows up in the value of H and e at horizon crossing, leading to a small and 
slowly oscillating tilt. The result reads (see also [27]): 



(CkCk') = (27r)35(3)(k + k') 



4M2je(t,)A;3 



a < 1 



(23) 



2.3 n-point correlator 

In order to evaluate the single vertex connected n-point correlation function we need to 
derive the n-th order interaction Lagrangian for vr. Since derivatives of the background carry 
potentially large factors of a (we will indeed see that the most interesting physics happens 
when a ^ 1), it is useful to make derivatives of the background quantities as explicit as 



possible by integrating by parts the action (11). We obtain 

1 



PI 



Tjln) ■n—l • 

-n^ 'It n 



(n-l) 



[n-l] 

Zj(n-—l)„n—l ■■ 



1^ 



TT 



1 



1)! 



6C 

671 



' ) 

osc/ 



where 



/n(vr) 



-1) 



2(n 



H 



-TT 



(24) 



(25) 



^In order to compare, notice that Eosc — ~Gbthere/ct- We prefer to use eosc rather than b because, as we 
stressed in the introduction, we find this notation to be more justified in terms of symmetries. Indeed, the 
constraint b < 1 that is often assumed in the literature comes from imposing the monotonicity of the potential 
for standard slow roll inflationary models. Such a constraint is not necessary. 
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and 5C/5Ti\^ is the linearized equation of motion. Here we have taken the leading order 
contribution in eosc, and used the following scaling for the time derivatives of H from (13): 

Hr^eH^ (26) 
a^^i/ ~ e eoscft""^^""^^ sin(a;t) for n>2 . 



The second term in the third line of (24) is suppressed with respect to the first by a factor 1/a. 
The terms in the last line proportional to the linear equation of motion of can be removed 
through the following field redefinition 



(27) 



This has the side effect of modifying the interaction Lagrangian £„ for all n > 4. However, 
the terms induced by the field redefinition are suppressed by at least one power of eosc ^ 1 
with respect to the leading interaction and so are negligible at the leading order. Nonetheless 
in conversion of the correlators of vr^ to those of vr these field redefinitions must be in principle 
taken into account as we will do later in this section. We are therefore left with the following 
leading interaction: 

Sn = - [ d^x ^M^, (//("+!) + SHH^^A < + ^(e^ ) . (28) 
J nl 

After canonical normalization of vr^ 

TT, = {-2Ml^Hf/^TTr , (29) 

the interaction Lagrangian becomes 

Cn{iic) = --Ml, (i7("+i) + 3MW) { —. ^ . (30) 

It is possible to show that this Lagrangian agrees with the n^^ order interaction Lagrangian 
obtained in |22] where a model of a slowly rolling scalar field was considered]^ The calculation 
is straightforward following [22] • The Fourier modes undergo resonance at k/a{t) ~ u and 
the integral can be done in the saddle point approximation. Using the conversion ( = —Hn, 



•^We will show the equivalence just at leading order in Cosc- The interaction Lagrangian of order n in [22] 
reads 

n\ 

Using the Friedmann equations that allows us to write F((/)) and (j) in terms of H and its derivatives 

V ^ M^i{H + 3H^) , 0^ = -2M^iH, (32) 



10 



we obtain 




i=l 



Bn{ki 



/ \n+l ^oscV27r 2n-7/2 / H' 

^ ' 4 \AtM, 



2 \ 



(37) 



2 



sin (a ln(ir/A;*)) + - cos (a \niKlK)) 

a ^ 



where K = ki. These expressions are vahd only for a ^ 1 and they agree with [22j. 



The resuh for a ~ 1 or smaller can be estimated by taking a — )■ 1 in (37). In this case 
we get an effective /nl ~ Cosc, which, after including the bound from the tilt of the 2-point 
function, results in a negligible amount of non-Gaussianity. We therefore safely restrict to 
the case a > 1 for the rest of the paper. We will explain in App. |A] and |B] the reason why 
coupling to gravity can be neglected and the linear conversion between ( and vr is sufficient. 



Let us first comment on the scaling of (37) with a. It is clear from (30) that as one 



goes to higher order interactions, at each level there is an extra time derivative which gives a 
factor of a; = aH. Moreover since n}^ order interaction Hamiltonian contains n fields with a 
wavefunction of the form 



H 



(1 + ikr)e 



{31 



in view of the fact that the resonance happens around k/a ^ u {or equivalently |A;r| ~ a), 
we expect to get two more factors of a at each order of {("')■ 

It is easy to get the n dependence of the formula above. We just need to realize that the 
ratio of (C")/(C^)"^^ scales as the ratio of the interacting Lagrangian Cn and the quadratic 



the derivatives of V with respect to 4> can be written in terms of time derivatives of H 



(33) 



(34) 



It is easy to see that in higher derivatives than V" , taking derivatives of (j) results in second or higher powers 
of Eosc- We obtain 



(35) 



and finahy using (32) and (311 



nl \ 



(-27fM2j)i/2 



(36) 



which is identical to the n order Lagrangian (30) for canonically normalized tt^ at leading order in eosc- 
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Lagrangian £2 evaluated at energy scale of order oj where the resonance happens: 



(e)"/^ ^ A ■ ^^^^ 

Notice the important difference here that the non-Gaussianities are dominated by interactions 
happening when the energy is of order uj contrary to the standard case when non-Gaussianities 
are dominated by interactions happening at energy of order H. This difference arises due to 



the resonance that happens when the mode is way inside the horizon. Plugging into (39) and 



using (26) we obtain 



' 'cjtt) ~ eosc [a Cj , (40) 



where we have used that vr at an energy scale E, tte, is related to vr at Hubble by he 



(E/H)7rH ~ E/H'^C and that C, ~ Htth. To get (37), there is another factor of a^/^ that 
comes from the saddle point approximation. More instructively, (39) can be also re- written 
as 



(C) 




n-2 



1/2 / \n-2 1/2 ( ^ 1 1/2 / "c 

Eosc « ' (wtt) ~ eosc a ' - — ■ , ,0 ~ Cose a ' 



n-2 



n-2 



~ eosca'^'[j;) ■ (41) 

where F is related to the unitarity bound Au of the theory due to the discrete-symmetry 
preserving terms by the relationship Au ~ AnF. Performing the same procedure keeping 
track of numerical factors, we find the value of F to be 

/ ■ \ 1/2 

cos (Uj(t + Tt)) = cos I t + -. TTc F ~ ^ ^ — . (42) 

^ ^ ^' ' (-2™2j)l/2 y iU ^ ' 

The perturbative series is defined as long as (C")/(C^)"'^^ ^ I5 which implies to <^ Au. 
This obviously corresponds to the case in which the resonance frequency is smaller than the 
unitarity bound of the theory: u ^ Au Thus the bound on the consistency of the effective 
theory implies that 

a < . (43) 

This is a very important bound that will have relevant observational consequences. Indeed 



eq. (40) is proportional to the signal to noise, and we see that within the validity of the 
effective theory the lowest n-order correlation functions have the leading signal to noise. 
Before being sure that this is the case, we need to take care of the mixing with gravity and of 



''Alternatively we could say that in the case oi ui ^ F, we could still perform the calculation of the n-point 
function but we should restrict ourself to frequencies much smaller than F. In this case the effect of the 
oscillatory functions would be suppressed by F/uj instead of leading to an enhancement. 
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the redefinition between C, and vr, and be sure that terms that we neglected do not represent an 
important contribution for higher order n-point functions. We prove that these contributions 
are neghgible respectively in App. |X] and [B] This allows us to conclude that as a general 
prediction of the effective field theory the leading signal to noise is in the 2-point function. 
It would be interesting to investigate if this conclusion can change if one considers specific 
UV completions of the Effective Theory where new degrees of freedom are included so that 
resonance frequencies larger than F can be consistently considered Fj 



3 Upper Bound and Observability 



The upper bound (43) from the consistency of the effective theory restricts the resonant 



n-point function amplitudes An in (37) to satisfy 



IA2I IA3I \An\ 



This implies that the resonant part of the 2-point function is the best observable prediction 
of this model, unless some (unexpected) degeneracies in the 2-point function are present 
This is one of the most important results of this paper. 

To derive this statement more rigourously, let us estimate the precision of a measurement 
of an n-point function using the Fisher matrix method (see for instance [28]). Suppose we 
measure Npi^ points in a volume V of the sky. Transforming to momentum space, it amounts 
to observing modes up to a k^ax momentum, where Vk^^^ ~ Npi^. The ratio of signal to 
noise can then be estimated as follows 

N^^^^^J ^ ^ y (2^---(2^ (Ck.---Ck„Ck.---CkJ ^''^ 

where for simplicity we have neglected the transfer functions and we have taken a 3-dimensional 
survey. Conclusions are not expected to depend qualitatively on these simplifications. Taking 



While a complete answer can only come from detailed studies of specific UV completions, we can anticipate 
that such a possibility seems unlikely to us. It seems indeed very hard to induce resonance effects at frequencies 
Lu > F, even after a UV completion has soften the amplitudes to preserve Unitarity at such energies, as in 
this case it is expectable that the oscillating components of the background should be erased by the vacuum 
quantum fluctuations of the fields when the modes have frequencies of order lo. On top of this, given that the 
role of the UV completion is to soften the amplitudes, it seems even harder the possibility that the eventual 
remaining signal would be stronger in the higher order n-point rather than in the 2-point function. 

^In [20], an analysis of the 2-point function was carried out for the CMB WMAP data, identifying a 
degeneracy with ilb, leading to weaker limits of at most a factor of 5 for some small region of values of a 
smaller than 10, when a specific value of the phase is chosen. These are much smaller values for a (by a factor 
of at least 10) than the maximum allowed ones. Since the signal in the 3-point function scales relative to one 
in the 2-point function as the square of the ratio of the a's, we conclude that this small degeneracy cannot be 
broken by analysis of higher order n-point functions. Further this degeneracy will most probably be broken 
completely by the Planck data or already at present by carrying out an analysis using the data from ACT 
and SPT where a higher number of acoustic peaks are measured with a high signal to noise. 
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(27r)^(5^(0) = V and neglecting for the moment logarithmic corrections, we obtain 



which after combining with (44) results in 



|(^(0)>|((C'»>--->|((C"»>---. (47) 

A careful computation, taking into account all numerical factors, confirms that the ratio of 
the signal to noise for the 3-point function compared to the one from the 2-point function 
reaches the highest value of 0.87 < 1 even when we push a to be so high to saturate the bound 



of (43) (see Appendix [C]) . For consistency of the perturbative expansion, a should be much 
smaller than asaturation, and the ratio of the signal from the 3-point function with respect to 
the one from the 2-point function gets suppressed by the ratio (a/ctsaturation)^ for smaller a's. 
Therefore the measurement of resonant part of the power spectrum is the most sensible test 
of these models. 

In the next section we will investigate generalizations of this resonant model by letting the 
additional operators in the Effective Theory that we neglected so far become relevant. We 
will try to see if they can lead to a significant level of non-Gaussianity with a larger signal to 
noise with respect to the 2-point function. The answer will be that this is not the case. 



4 Generalization 

Within the EFT of Inflation, one can easily generalize resonant models to incorporate small 
speed of sound or large inflaton self-interactions. It is enough to add higher order geometric 



operators to (11), which after expansion in terms of vr (restoring gauge invariance) result in 
new interactions. The coupling coefficients in these new interactions will naturally have small 
oscillating components which generate resonant non-Gaussianities similar to ones studied 
before. However, there is an additional way in which resonant effects can be important. So 
far we have simply considered the effect of an oscillating coupling on a standard Bunch- 
Davies wavefunction as the coupling was already of order eosc- Here this correction to the 
Bunch-Davies wavefunction will need to be considered as the interactions are not necessarily 
suppressed by eosc anymore. 

We will show in the following that the conclusion of the last section is generically true 
in these two additional cases, namely, the higher order correlators remain smaller than the 
modificatons of the 2-point function as long as the effective theory is valid and natural. 
However there is an extra enhancement of the bi-spectrum in the large folded limit in the 
second case above. As we will see, this will make the signal-to-noise originating from non- 
Gaussianities larger than the one in the 2-point function only in a very small marginal region 
of parameter space where the bound on a is saturated. 
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4.1 Oscillating Couplings 

In this subsection we show that the oscillating coefficients of self-interactions of tt (if present) 
are tied by quantum loops to the modification of the power spectrum in such a way that the 
contribution to the signal to noise from higher correlators is always suppressed, at least if we 
do not consider fine-tuned theories. 

Let us consider the simplest scenario as a benchmark. Substituting (16) into (18) the 
linear equation of motion for vr schematically takes the form 

TT + 3H{1 + Cosca cos(wt)) tV - 9V = . (48) 

This modification of the equation of motion excites a negative frequency part in the mode- 
functions [201 [21] (this will also be discussed in more details in ^4.2.1) which leads to 

5(C^>~e_«^/^(C^). (49) 

Now consider adding the following cubic interaction 

(l + Acos(a;t))^, (50) 

where A is the high energy cutoff of the theory at which the interaction tt^ becomes strongly 
coupled. If we let loop corrections run until the strong coupling scale A (a necessary condition 
if we do not give an explicit UV completion cutting off the loops and making them convergent), 
the contribution of this coupling to the 1-loop renormalization of the kinetic term is of order 

(l + Acos(u;t))^7T^ (51) 

Neglecting numerical factors, now the linear equation of motion becomes 

7f + 3/7 (1 + (cosc + A)a sin(wt)) ir - d\ = . (52) 

We see that the loop correction effectively renormalized eosc- In order for this renormalization 
to be at most of order one, so that the theory is technically natural, we need to impose 

A ^ eosc ■ (53) 



The level of resonant non-Gaussianity produced by the second term of ( 50 ) is easy to estimate 
following sec. 2.3 We obtain 



(^2)3/2 ^^osca^' (Xj ' ^^^^ 

where we have inserted a factor of a^^^ coming from the saddle point approximation. This is 
always smaller than (49) within the range of validity of the effective theory w ^ A. In the 
computation of the non-Gaussianity at leading order in eosc, we should consistently also include 
the contribution from the perturbed wavefunction applied to a non-oscillatory coupling. We 
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will describe this contribution (which turns out to be the leading one) in the next section. 
Here we concentrate only on the non-Gaussianity resulting from considering the oscillatory 
coupling on the unperturbed wavefunctions. 

The arguments above can be easily generalized to any other interaction. For instance 
the resonant 4-point correlator in the presence of a large quartic interaction tt^ as in [5] is 
constrained to 



((2)2 ^ VA 



^eosc«'/'(T) . (55) 



We will complete the study of this case in sec. 4.2.4 Thus, we conclude that by adding 



an oscillating part to the couplings one cannot escape the observability bound of previous 
section. 

Notice that in doing these estimates we have always used the cutoff of the theory and the 
canonically normalized tt. By following the same steps as in [5], we can easily see that even in 
the case of small speed of sound we do not obtain an enhancement of the non-Gaussianities 
because for any Cg the non-gravitational sector of the Lagrangian can be made effectively 
Lorentz invariant by rescaling and defining tTc = {—2M'^yHcsY^'^ti and reabsoring the 

remaining factors of Cg in a rescaled cutoff A 0. Notice that now there is a stronger upper 
bound on a coming from imposing u to be within the regime of validity of the effective theory. 
In this case the UV cutoff is given by [1] 

A^ ~ IQtx'^ Ml^Hcl , (56) 

which after using (C^) = H^/AMp^Hcg yields 

A / 2n 



'^saturation (Cs) Cg ^^^2^1/2y ~ '-'s'^^saturation ('^s 1) • (^''') 

Since the ratio of the signal to noise of the 3-point function with respect to the 2-point function 
goes proportionally a^/c^, we explicitly see that taking small Cg does not help to increase the 
relative signal to noise ratio for the 3-point function Q 

This conclusion can also be extended to multi-field inflation. In the case when some 
additional light fields (a) have an effect on the duration of inflation, there are two possible 
sources of large non-Gaussianity [9]. The first is due to their self- interactions. In this case 
our analysis of single field inflation applies equally well to correlators of a's. However even 
with negligible self-interactions, non-Gaussianity may be generated because a fields might 
affect the duration of inflation in a non-linear way. This is the second way in which non- 
Gaussianities can be generated: that is through the coefficient of proportionality between ( 
and (j2. If d'^C^jda'^ is oscillatory, one might hope to have oscillatory non-Gaussianities of 
the local kind. However, it is very unclear how one could imagine not having a comparable- 
size oscillatory component in the linear coefficient dC,/da, which would lead to an oscillatory 
component in the 2-point function that dominates the signal. 



^Later we will see that there is an enhancement of the relative signal-to-noise coming from the folded limit 
of the triangles by a factor of order ^a. So even in this case the maximum relative signal to noise decreases 
by reducing Cj. 
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4.2 Resonance from Corrections to the Mode-functions 
Shapes 



Folded 



Another source of resonant non-gaussianity in theories with large deviation from = 1, or 
more generally in theories with enhanced cubic derivative self-interactions, is the correction 



to the Bunch-Davies wavefunction due to the oscillating background (see e.g. (48)). This 
correction basically amounts to the addition of a small negative-frequency part to the wave- 
function after the momentum matches the resonance frequency. In the original setup without 
large derivative self-interactions, these corrections could be neglected because the cubic and 
higher order interactions generated by the expansion of cos uj{t + tt) are already proportional 
to the small amplitude of oscillations eosc, and so this effect would be of order e^g^, ^ eosc- 

We will show now that in the presence of large interactions these contributions to the 
non-Gaussianity have a comparable size with the previously considered cases (i.e. {(^) ~ 
Cosc^^'^^C'^)) aiid thus they should be taken into account. Even more importantly, it is further 
enhanced in the limit of folded triangles. This is a general feature of non-Gaussianities 
generated through the modification of the vacuum wavefunction [29], which have already 
been calculated for resonant models in by means of approximating the coefficient of the 
negative frequency mode by a smoothed step function. Here we avoid this approximation and 
rely instead on the saddle point method. We recover the same scaling, namely (C^) oc a^^^ 
away from the folded limit and (C^) oc a^/^ in the folded limit. We emphasize that such an 
enhanced 3-point function signal can potentially escape the observability constraint of the 
last section when a is close to its upper bound. 

Following [21], we will first review the derivation of the modification to the Bunch-Davies 
vacuum and then compute its contribution to the resonant non-Gaussianity. Later we present 
the results of numerical computation of the signal to noise ratio, showing that despite the 
enhancement in the folded limit the 3-point function remains subdominant in practically the 
whole relevant region of parameter space. The reader who is not interested in the technical 



details can refer directly to ( 78 ) and subsequent equations for final results. 



4.2.1 Correction to the Wavefunction 

Let us consider the effect of adding the following operators 



2! 



(5^ 



to the minimal action (|11|). After restoring diff. invariance the quadratic action for vr becomes 



S 



(59) 



where 



Ml^H + 2M| 



(60) 
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The linear equation of motion reads 

- ^^^^^.VTk + TTk = , (61) 

where we defined x = —Cgkr and 6 is now generahzed to 

c? d { H 



2HH dt \ c 



2 



(62) 



This equation is identical to the one studied before in ( 19 ) . We assume that 6 has an oscillating 
component 

Sosc = - ^eoscO sin ut , (63) 

as used to be the case when Cg = 1, and it is a generalization of the definition of eosc to the 
Cs 7^ 1 case. The constant part of 6 can be safely neglected as the resulting effect is as usual 
slow-roll suppressed. 

This oscillatory component excites the Bunch-Davies vacuum [20], and the change in the 
usual vacuum solution to first order in eosc can be written as a small admixture with negative 
frequency mode 

7rk(x) = vrf K(a;) + 4"^(x)iz_(x)] = Trf [(1 - 2x)e"' + 4"^(x)(l + tx)e-'^] , (64) 

with TT^^ = (y/4MpjeCsA;^/^)~^/^ . Substituting (64) in (61) and solving perturbatively in 5, 
we obtain at linear level El 



c 



^ ^ ' n + ix'V J x"^ ^ ' 



Clearly 7rk(x) freezes at late times which is evident from the quickly vanishing dc\. ^ jdx ^ x 
for small x. The expression above can be simplified for the resonance period which takes 
place well inside the horizon x ^ 1. Using (63) and noting that lot = — aln(x/A;) plus a 
/c- independent constant, one obtains 

ci'\x) ~ ieosc / dx'e'^'''' cos{ut{x')) . (66) 

J oo 

As thoroughly discussed in |20j, c^i^\x) starts from zero at early times {x — j- oo), jumps at 
resonance period x ~ a/2 and stays constant afterwards. Both the step-like behavior and 

^The explicit analytic solution can be found and the correction to the unperturbed wavefunction reads: 

Stt = -^^^^^—^e~"-"i-,^Hr- ({a + zf + ae^^'^^k^ + i) ({a - zf {-i^Hr^ E.^W^zkr^) 
^JAMl^eCs 4 [a^ + 1) k'^'^ V V 

- (a + if E.,a+i{2ikr])^ + [a - i){--qHf'°'{a - iak-q + krj - i) ~ i (a^ + l) fcr/^ , 
where E^(x) represents the exponential integral function of index v. 
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the small oscillations on top of that are important to obtain the resonant behavior and the 
additional enhancement in the folded limit. 



Expression (66) gives a good approximation for x not much smaller than one and so 
it covers the resonance period x > 1. Therefore we will use this simplified expression in 



the calculation below. Let us note that (66) will make some of the terms, for example the 
contribution to the 3-point function from the vertex 7r(c}i7r)^, IR divergent. This is a fake 



divergence which will disappear once one uses the full expression (65) which is different from 



(66) for small x. Equivalently, to get the correct result at leading order one can simply 



use (66) throughout and simply disregard the contribution from the small x region. 



The value of cl ■* at late times, when x — > can still be well approximated with the saddle 



point approximation for (66), since the phase is stationary at x ^ 1. We obtain 



4 



and the power spectrum can at this point be easily derived to obtain (20) 



(67) 



4.2.2 3-point function 



After expanding ( 58 ) in terms of tt the cubic action reads as follows 



{61 



with X^sMp^Hc- 



' ^ -MlMcf - 1) - -^Mt , 



where as discussed in [T] M3 is a free parameter. It is technically natural to choose M| ~ 
M^/c^g ~ —MpyH/c'^, so that tt^ gives the same level of non-Gaussianity as 7r(9j7r)^. This 
choice corresponds to X^^s ~ 1. As mentioned earlier, in the presence of large cubic self- 
interactions, modification of the wavefunction can generate large non-Gaussianities. At lowest 
order it suffices to use the negative frequency term, instead of the positive-frequency one, for 
one of the fields in the calculation of the 3-point function using first order perturbation theory. 
For instance from tt^ we get 



{-i)Xi,3MlyHc- 



n 



dr 



{—)* * \/ * / * / 



-|- perm. -|- c.c. -|- other choices of c*- ^ 



(69) 



where m+ and u_ are defined in (64) as the positive and negative frequency solutions and 



we have dropped their arguments to avoid notational clutter. Prime denotes derivative with 
respect to the conformal time r. It can be straightforwardly checked that substituting the 
negative modes in the final wavefunctions gives a negligible effect because there is no resonance 
enhancement. 



Because of the oscillating behavior of c 



(-), 



X 



one discussed in the previous subsection. This explains the scaling (C^) ~ eosctt^'^^(C^)^ away 



see ( 66 ) ) the integral is very similar to the 

3\ 
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from the folded limit. The further enhancement in the folded limit is due to the fact that in 
this limit the phase of the integral is almost constant and c_ can be taken to be a step function 
that rises from to eoscCt^''^ at resonance frequency kr ~ —a. This effectively enhances the 
time integral by a factor of a. 

Now we will show that since the resonance happens well inside the horizon, when |A;r| ^ 1, 
in the integral 



° dr 



T 



(—)* * \/ * f * / 



we can keep only the terms with highest power of r. Let us use the variables 

Xi = CgkiT , and yi = \ — — ^ ^ kj , 

which are useful because y — >■ 0~ corresponds to the folded limit. Ij^-i can be written as 



(70) 



(71) 



where we introduced 



r3 1.2 1.2 

ki 



(/i(2)(yi) + /n«(yi)-^/ii^'nyi)) 



(72) 



rfxix^e'"i^^c(-)*(xi) 



dxi 



(73) 



We only know an integral representation for K^i) and therefore it is convenient to integrate 
by parts to obtain 



-I 



dvi [ iyi 



m Jo 



dxi(e 



txiyi 



(74) 



Since c^"-* is zero at early times, c^^\x — )■ oo) = 0, the first term in (74) vanishes. With help 
of the explicit form of dc^^^ /dx which can be extracted from (66) the remaining integrals can 
be calculated using the saddle point approximation. 

The leading order in a contribution comes from the highest derivative which acts on yi in 
the argument of the exponent: 



1/2 / — 
. eoscft yTT ia ln(fci/A:*) 

t , — C- 



V2 

and similarly 



j{n) _ ^(n) _ 
-'II ~ -^11 — 



Qn g-ialn(2-j/i) ^rig-io ln(2-j/i) 

dy^ (2-yi) - (2-i/i)"+i 



(75) 



-I -1 



-ia ln(ki/kt) 



V2 



j{n) _ ^{n) 



Qri lQ-ioi\n(2-y) 



dyi yi{2-y) 



y=yi 



%a e 



rip-ia ln(2-j/i) 



l/i(2-|/i)«+i 



(76) 
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The approximation (75) is valid for any yi while the last passage of (76) is only valid for 



sufficiently large \yi\ > 2na ^ far away from the folded limit. In the folded limit ?/i — ?■ we 
can use the expression for X^"'^ as it approaches yi — t- 0~ 



(n) 



Q,n+lg-io In 2 



— ^ 



1 + 



iayi{n+l) a'^y1(n+l) 



2{n+2) 



8(n+3) 



+ . . . 



(77) 



2"+2(n+l) 

In Figure [l] we compare the explicit expression for Xp-* with the approximations ( 76||77 ). For 
large a the approximate expressions give a good agreement for both yi — )■ 0~ and yi < —1/a 
regions. In the region yi ~ —1/a the exact expression for Xj"'' oscillates with a rapidly 
changing amplitude and both approximations break down |^ 




(2) 

Figure 1: Explicit expression for Xj for a 
and small yi correspondingly. 



100 and the approximations (76 77) for large 



Since Jj"] ~ a" we conclude that only X*^") with the highest n are important in (72), 
confirming that we can keep only the terms with the highest powers of r. Substituting the 

(2) (2) 

expressions for Jj and Jjj , summing over permutations and using ( ~ —Htt we finally find 
the correlator of three C's away from the folded limit 

<C^>,3 — 3y/|A.3C-A?eo..a^/^ ^ " ^) k^kk^^ sin(alnir/A:.) , (78) 

with K = ki + k2 + ks and we have introduced = / (AeCsMp^) . In the folded limit 
ki ^ k2 + k^, we have 

{('YT' ^ ^Xn^c;'AlV2^eosca'/'k%k, cos{aln2k,/k^) . (79) 

As anticipated, we see that away from the folded limit we have the same scaling as before 
~ a^/^, while at the folded limit we have ~ a^/^. This expression holds approximately for a 
region in yi ~ [— l/a,0]. 



^This means that for yi —1/a, which corresponds to a quasi folded hmit, we should use the exact 
expression, but at this point the formulae become bulky and we do not compute them explicitly. However 
from the plot we see that this region corresponds to simple interpolation between the two extreme regions 
where the expressions are simple, without any important change. 
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Figure 2: The shape of non-gaussian signal (78) induced by vertex tt^ plotted for a = 50 



away from the folded limit as a function of Xi = ki/k\. 



The computation of the contribution to the 3-point function due to the vertex TcldiTT)"^ is 
similar 



■3 1 ' 



-1 



I' + perm. + c.c. + other choices of c*- ^ 

(80) 



I' = Cs 



fc2fc3(k2 ■ kg) (2) (2). 

-[II ' [yi) - ill^' [yi)) 



-C.S^^/f)(.2)-C. 



A;?A;2(k2-k3) 



k1 



^1) 



Here we dropped the subleading Jj jj with n < 2. Now it is straightforward to get the 
contribution to the 3-point function 



(a,7r)2 



7^(Cs - Ij^^^eosca p 7,2 7,2 /^-S 

Z( Ay-j^ '1^2 V 



X 



^2) 



^ ^IZ^i'-I]^*^!) f 1 + 5^- 1 1 sin(alni^/A;,) 



away from the folded limit, and 



- :^(c- - l)A?x^e„.aV^^eos(aln2fc,/fc.) 



r(9«7r)2 32 

in the folded limit fci — )■ ^2 + /^a 



^3) 
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Figure 3: The shape of non-gaussian signal for the 3-point function generated by the vertex 



7r(9j7r)^ (82) plotted for a = 100 away from the folded limit. 



For completeness, we also give the shape induced by the oscillatory couplings on the 
unperturbed wavefunction. This is given by [SU] 



osc'-* 



sin(Q; In K/k^) , 

(84) 



where A^3,A^(g.^)2 ~ 1, correspond to tt"^ and n{diny interactions, respectively. As described 
in the former subsection, we see that this contribution is also subleading in the sense of signal 
to noise ratio with respect to the 2-point function. 

4.2.3 Upper Bound and Observability for the Folded Shapes 

The shapes of the 3-point functions we computed are very different from the standard ones 
that are analyzed. Indeed, because of the oscillations in fc-space, in App. [D]we show that for 
reasonable high a's they become orthogonal to the standard shapes. This means that current 
bounds on non-Gaussianities do not relevantly constrain these shapes. However, we are now 
going to see that the leading signal-to-noise for these models is in the 2-point function. 

Because of the enhancement in the folded limit, the signal-to-noise ratio can be enhanced 
with respect to the naive one. Indeed the naive inequality {(^) /(C^)^''^ < 5 (C^) /(C^); is not 
satisfied for configurations close to the folded limit, and so one might wonder if the enhanced 
signal-to-noise in the 3-point function coming from those configurations can potentially make 
the 3-point function the leading observable. The shape in the folded limit was given by 



eq. (78). This shape is valid for a region of parameters space from ?/ ~ 0~ to y ~ 1/a. The 
3-point function is enhanced in this region by a factor of a so that the signal-to-noise gets 
enhanced just proportionally to ^/a. 
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Figure 4: Ratio of the signal to noise ratio for the 3-point function generated by vertex tt 



^3 



(78, 79) to the signal to noise ratio of the resonant correction to the 2-point function (20) 
plotted as a function of a for = 0.5 on the left, and = 0.1 on the right. The shaded 
region corresponds to the values of a where the effective theory breaks down, and the allowed 
values of a should be well within the white region. We see that the region where the signal 
to noise from the 3-point function is bigger than the signal to noise from the 2-point function 
is irrelevantly small and corresponds to when the theory is not under parametric control. 



In Fig. |4] we plot the ratio of the signal to noise for the 3-point function versus the one 
for the 2-point function. We see that there is a very small region in which this ratio is larger 
than one. This happens only when a becomes extremely close to its upper bound. Notice 
that in order to have a controlled perturbative expansion the parameter a needs to be much 
smaller than that. We therefore conclude that the leading signal to noise comes from the 
2-point function. 



4.2.4 4-point function 

In [5], it was shown that there are technically natural models of single clock inflation where 
the 4-point function has a larger signal to noise than the 3-point function. For models 
not very close to de Sitter space, the 4-point function is induced by only one operator 
M| (1 + /io cos(a;)f:)) TT^. Following the steps of [5J, we can rescale the spatial coordinates 
and canonically normalize vr to show that this term takes the form 

71-4 

(i^xdta^ (1 + /io COS (wt)) , (85) 



where tTc = (— 2ijM|iCs)"^/^7r, x = x/cs, and 



(hm, 



2 



- ' ..,7 ■ (86) 



Mfcs 
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Let us impose that the signal induce by the 4-point function is bigger than the one from the 
3-point function. This is achieved by imposing 



> 1 



UJ^Cs 



3/2 



^7) 



In [5] it is shown that it is technically natural for Cg to be order one or smaller. The result 
above further implies that the maximum a that we can have is (compare to ([57))) 



A 



'^^saturation 



UA 



H 



'-^saturation 



(^2)1/4 • 



If we now turn to coefficient of the oscillating term fiQ, it should not renormalize the oscillating 
kinetic term by a large amount, implying the constraint 



^9) 



The leading contribution to the four point function consists of two parts. First, there is a 
contribution due to the correction to the Bunch-Davies wave-function which is fio independent 



11 AM^,ec,kf 



7 I ( — )* * \/ * / * / * / 



+ perm. + c.c. 



(90) 



Using the variables 



Xi — —CgkiT , and |/j — 1 — — ^ ^ kj 



the contribution above can be rewritten as 

4 



n 4A/f2,, 



i=l 



AM^^ec.kf 



(91) 



kf 

+ perm. + c.c. , (92) 

where we neglected the subleading contribution of Jjjj with n < 4. With the help of (75) 
and (76) and using ( —Hn, we obtain 



(93) 



away from the folded limit, and 



4, \ folded 



3A/2^en.ra^^/^ 1 



160 klk2k^ki 



cos(q; In 2ki/k-^) 



(94) 
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in the folded limit ki = k2 + + k^. We see that there is an enhancement proportional to a. 

Second, there is a contribution due to the oscillating coupling /xq cos(a;t) which is also easy 
to calculate 



\ / ^toTT' 



n 





'OO 

+ C.C. 



7 */ */ */ / x\ 



(95) 



After the usual change of variables x = —CgKr this reduces to the integral which can be 
calculated using the saddle point approximation 



\^ / itnvT 



Ml 



1 



sin (a \nK/ki,) 



(96) 



The explicit calculation confirms the qualitative conclusions of section 2.3 the four point 



function scales as a^^^A^. In the folded limit there is an enhancement of oder a. 

Let us consider the ratio of the signal-to-noise with respect to the 2-point function. If we 
ignore the contribution from the folded limit, we obtain 



Ml 



A 4 ~ 



< 1 



(97) 



Therefore the signal to noise ratio of the four point function will be small. The folded limit can 
give a further enhancement that is proportional to a^/^ as for the 3-point function discussed 
in the previous section. We expect the numerical coefficient for this scaling, that was very 
small for the 3-point function, to be even smaller in the case of the 4-point function because 
we expect that the fraction of parameter space affected by the folded limit to be numerically 
smaller in the case of the 4-point function. Because a^^^ is a slowly growing function we 
expect the value of signal to noise to become large (compared with the one for the 2-point 
function) only in the marginal area of a ~ asaturation where the effective theory is becoming 
strongly coupled. 

We finally add a comment on the multifield infiationary case. As discussed in in 
the case of multifield infiation it is possible to have a quartic operator of the form (Sjcr)^, 
with a being an additional scalar field, without the presence at the same time of any other 
cubic or quartic terms. This is impossible in the context of single clock infiation |BJ. If we 
suppose that the coefficient of this operator or of the quadratic Lagrangian have an oscillatory 
component, then we may ask ourselves if this quartic operator can generate large oscillating 
non-Gaussianities. The answer is again no. Indeed, under renormalization, this operator 
would induce an oscillating component in the part of the quadratic Lagrangian proportional 
to (Sjcr)^, without renormalizing the time kinetic part. The resulting linear equation of motion 



would be different than the one we studied in sec. 4.2.1 with the oscillating term simply sitting 
in front of the spatial kinetic term. Apart for this detailed difference, the effect of the 2-point 
function is of the same order as the one found before, as it can be verified by explicitly finding 
the solution for the wavefunction. So, we conclude that also in this case the leading signal is 
on the 2-point function. 
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5 Conclusions 



Non-Gaussianities provide one of the most important probes of Physics of the Early Universe. 
They are closely connected to the interacting structure of the inflationary Lagrangian and 
their signal is not strongly constrained by the symmetries of the problem, as it happens for the 
2-point function. This leaves room for a large amount of non-trivial information to be encoded 
in higher order correlation functions. For these reasons their detection would represent a huge 
improvement in our capabilities to probe inflation. It is therefore of the utmost importance 
to understand all possible non-Gaussian signatures that can be generated during Inflation. In 
this context the Effective Field Theory (EFT) of Inflation represents the ideal setup, as it is 
particularly well-suited for exploring the whole parameter space of Inflationary models, and 
classifying them in terms of symmetries of the underlying Lagrangian. The EFT of Inflation 
also naturally provides an insight into the scaling and the relative importance of the various 
operators. Different theories can indeed be classified according to the different symmetries 
they respect, very similarly to what happens in the particle physics context. 

In this paper wc have considered the case in which the Goldstone boson of time trans- 
lations is endowed with a softly broken discrete shift symmetry, and we have studied the 
observational consequences of this pattern of symmetry breaking. Our study unifies and gen- 
eralizes previously known infiationary models such as Axion Monodromy Infiation. Because 
of the non-linear realization of time-diff.s, the discrete shift symmetry forces the presence of 
parameters in the EFT Lagrangian that oscillate with time. This induces a resonance as the 
modes expand inside the horizon. The resulting leading signature is the presence of oscillatory 
features in the power spectrum of density fiuctuations. 

Additionally, the same symmetries force the presence of interacting operators for the 
Goldstone boson that, when combined with the resonant effects, lead to non-Gaussianities 
with remarkable oscillatory features. Within the regime of validity of the effective theory and 
for technically natural theories, we find that the signal to noise ratio of this non-Gaussian 
features (which represents how large these non-Gaussianities are) is always smaller than the 
one in the 2-point function. That is to say that analysis in search of these oscillatory features 
can be safely restricted to the 2-point function. Higher-order correlation functions induced by 
oscillatory operators can be detected only if the oscillatory 2-point function is detected with 
high signal to noise, or if some (unexpected) degeneracies are present in the 2-point function. 
Such a detection would be a remarkable confirmation that Infiation actually took place, and 
that the Goldstone boson was endowed with a discrete shift symmetry. 

Our bound originates from the fact that within the Effective Theory there is a maximum 
energy scale h.u beyond which unitarity is violated. Non-Gaussianities are small because 
resonance frequencies larger than Kjj cannot be considered. It would be interesting to see if 
our conclusions remain unaltered when one considers specific UV completions of the Effective 
Theory, where, probably after introducing new degrees of freedom, one could explore regimes 
of energies larger than Ku. As we briefly argue in the main part of the paper, the possibility 
that our conclusions might change seems unhkely to us. A definitive answer can however 
come only from studies of detailed UV completions. 
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A Mixing with Gravity 

In all the calculations we have always concentrated on the decoupling limit and we neglected 
the mixing with gravity. This was due to the intuitive fact that the resonance frequency 
was way inside the horizon when non-Gaussianities are large: a ^ 1, and so we expect that 
mixing with gravity should be suppressed either by the ratio of these two scales, which is 
indeed 1/a, or by the slow roll coefficients. In order to check our intuition, we write the 
metric using the ADM parametrization 

ds"^ = -N^dt^ + Qij [dx' + N'dt) [dx^ + N^dt) , (98) 

we neglect tensor modes, and we choose the spatially flat gauge 

Qij = a^5ij . (99) 

For simplicity we neglect tensor fluctuations because they will be negligible. Now we should 
solve the constraint ADM parameters iV* at all orders in vr and plug them into the first 
line of the action ([T]). In this gauge the action reads 

S = j d'xa^Ml,i^-^{El-E') (100) 

+ H{t + tt) [-N-'^{1 + iif + 2N-\l + 'k)N'diTi + N {diTif - N'^ {N' d-ijif] 
-N{'iH^{t + iT) + H{t + u))^ , 

where 

Eij = - = a^HSij - N^ij) , N^^.j) = ^{N,.j + A^,.,) , (101) 

and ; stands for the covariant derivative with respect to the spatial metric Qij, which in this 
gauge is simply an ordinary derivative. In this gauge, the constraint equations read 

dj[N-\Ei - 6iE)] + 2N-^H{t + 7r)[(l + 7r)9i7r - N^djTrdin] = , (102) 

- ^(4 - E') - {3H'{t + vr) + H{t + tt)) 
+H{t + 7r)[Ar-2(l + 7V)2 - 2N-\l + ■k)N'di'K + {diTif + N-'^^N'Omf] = . (103) 

10 We remind that to obtain the action for tt at the n*'* order we just need to solve for the N, Ni up to 
order (n- 2)|3T]. 
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We are interested in solving these equations perturbatively in vr. It is therefore useful to 
re-arrange the equations in a perturbative fashion as follows 

SN^i = -N-^^^[{l + n)d,Ti-N^d,Tid,Ti] 

AA2 . 1 . . 

dAN-\5lNl--im + Nf)], (104) 



2H(t) 

where we have defined = 1 + 6N. Note that only the terms that come from the the time- 
dependent coefficients of tadpoles y/—gSg^^ and y/—g have an explicit t + n dependence as 
the Einstein-Hilbert action is invariant on its own. The first order solution is 

SM^) = -fvr , (1), = ^dt{m) . (106) 



It is not easy to solve equations (104) and (105) beyond linear order. However it is straight- 
forward to understand the order in a at which they contribute. From the structure of the 
equations we see that at each new order we get only one extra derivative. This extra deriva- 
tive either comes from explicit derivatives of vr, or from Taylor expansion of H{t + vr) and 
H'^{t + vr), or substitution of lower order solutions for N and iV* in the r.h.s. The leading 
^th order solution therefore scales as 

5Ar(") ^ («) ^ eeosca""^^"vr" . (107) 

Interaction terms due to mixing with gravity will originate from substituting the solutions 
for N,N^ into the action (100). At a given order in vr, we can compare the interactions 
due to mixing with gravity with the leading one of order M|ji7*^""'"^W"'. Going through all 
interactions is quite tedious, and we simply consider the two leading terms. One comes from 
the interaction operator associated with 



£(9,^^)2 ~ Ml,{d,N,Y ~ Mil 7r"-2(9,7r)2 , (108) 



^(m)^(n-m) 

which is suppressed with respect to the leading one 

~ M2ji7("+i)7r" , (109) 

by the factor 

^^!^~»„.<.^ (110) 

Another interaction that scales parametrically in a different way comes from terms of the 
form 

^Hit+.)SN-^HH(-h- , (111) 
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where we have taken the hnear term in 6N as the non-hnear terms in 6N would give rise to 
further suppressed interactions. This compares to the leading term as 



Finally, it is possible to check that diagrams where the graviton is exchanged, are even more 
suppressed than the one we consider. We therefore conclude that mixing with gravity is 
suppressed by Min(e/a, e eosc) ^ e- 



B TT to Conversion and Field Redefinition 

In the EFT of Inflation we computed n-point functions of the Goldstone boson vr, which 
is the field that naturally manifests the decoupling limit. However, we are interested in 
correlation functions of the curvature perturbation C,, which has the useful property of being 
time-independent for modes outside the horizon. In order to convert vr n-point functions into 
C n-point functions, we need to find the non-linear relations between tt and C, at the required 
level. In this Appendix we show that only the linear relationship C, = —Hn is important for 
n-point functions with detectable signal to noise ratio. 



B.l The non- linear relation 

In order to obtain the non-linear relation between tt and (, we need to perform a gauge 
transformation from any gauge where vr is defined to the unitary gauge where ( is defined 
and where tt is set to zero. Notice that apart for terms coming from the mixing with gravity 
(which are subleading), we could choose to describe tt in any gauge distant enough from the 
unitary gauge PQ |32] • Explicit results will differ by subleading terms coming from the mixing 
with gravity. We can choose for example to define vr in the spatially flat gauge: 

gi, = a{tf5ij . (113) 

Here and throughout we neglect tensor modes, as they do not change the answer at leading 
order. In the unitary gauge tt = and the scalar degree of freedom is encoded in the metric: 

= a{tfe^^5,, . (114) 

In order to set vr = 0, we perform a time diffeomorphism 

i = t + T{t,x), (115) 

where t is the time variable of the vr-gauge, t is that of the unitary gauge and T(t, x) can be 
found iteratively from 

7r(t,x) +T(t,a;) = 0. (116) 
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For example this gives at second order 



T = -TT + TTTT + C(7r^ 



(117) 



In the ADM parametrization, the metric in the spatially-flat gauge is of the form 
ds'^ = -N^dP + a^{i)6ij [dx' + N'dt) [dx^ + N^dt) , 



:ii8) 



where N and A^* are given by equations ( 104[l05 ). Applying the diffeomorphism (115) 
to (118), the spatial part of the metric becomes 



Qij = o?{t + T) 



S,. 



+ diTdjT 



1 



a2(t + T) 



(119) 



where Ni = N^. This metric is not of the form of ( 114[ ), due to some second and higher 
order terms in tt. Therefore one is forced to make an additional spatial diffeomorphism 

X* = X* + e^{t, x). 

It is very difficult to derive explicitly the full non-linear spatial diff. necessary to reach 
the C-gauge. However, since we need to convert vr fluctuations to ( fluctuations only when 
all modes are well outside the horizon, we can restrict ourselves only to leading order in the 
spatial derivatives. This simplifies greatly the derivation because, as we can see, the last 
extra terms in (119) vanish in the limit that the spatial derivatives for all fields vanish 



^^While the terms in T do not contain any spatial derivative, this is not the case for Ni, which must start at 
least at first order in spatial derivatives. We are going to prove by induction in the order of the fluctuations 
this quite intuitive fact. We can imagine to perform the diff. ai* = a:* + e^{t,x) and expand perturbatively 
in fluctuations: e] being of order n in the fluctuations. Notice that starts at second order. Then we can 



An} 



5k 



Indeed, similarly to what is done in [c 



], the condition that e|"^ needs to 



imagine to solve iteratively in 
satisfy is 

' ' ^ (120) 

where 6h\j^ represents the n-th order part of the last three terms in (119) and where . . . stands for terms 

involving lower order terms in : e|" , e|" ... or analogously in /3. By taking the trace and by 
applying didj to (120), we can solve for to be schematically 



/3 



{"} 



6h 



did. 



Q2 ""-^S 



(121) 



where . . . contains lower order terms and where we have neglected numerical factors since it will turn out 
that we are not interested in the actual explicit solution. Expanding in smallncss of the spatial derivatives, 
if we assume that the lower order terms have the property that TV* and are at least of first order in the 
spatial derivatives and /? is at least of second order, then we s ee th at the n*'' order /3^"^ is of second order in 



spatial derivatives as well, and similarly upon substitution in (|120| we conclude that also e\"^ is of first order 



(n) 

in spatial derivatives. With this information we are now going to show that also start at first odcr in 

the derivatives. Indeed in [31] it is shown to all orders that A^f'''' in C gauge starts at first oder in the spatial 
derivatives, and upon performing the spatial diff. parametrized by e^, it is easy to see that at this point also 
Ni in TT-gauge starts at first order. From this we see that assuming that at the orders 1, . . . , n — 1, TV^ and 
are of first order and P is of second order in spatial derivatives, we obtain that the same condition holds for 
71*'' order. Since in ( 106 1 we show explicitly that this counting of derivatives holds at the lowest orders in 



the 
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This implies that in this hmit the spatial diff. is irrelevant, and the relation between vr and 
C reduces to 



which can be combined with (116) to obtain, for example up to quadratic order 

( = -Hn + ^Hn^ + Hnic + ^(Tr^) . (124) 

Notice that due to the mixing with gravity vr is time dependent even after horizon crossing. 
By using the relation ( = —Htc and the fact that ( is constant outside of the horizon, we 
derive tt ~ eHn. This tells us that the time dependence can be neglected at leading order 
in the slow roll parameters. We conclude that for our purposes the effect of the non-linear 
relation between ( and vr has the following scaling 

^ if{"-i)vr« ~eeosca""^C" , (125) 

where C^"^ stays for the contribution of order n in powers of (. 



B.2 Corrections due to non-linear tt to ( conversion 

We are now going to show that the effect of the field redefinition arising from the oscillatory 
part beyond linear order is irrelevant for n-point functions that have a non-negligible signal 
to noise. In order to do this, we can use the constraint eosc ^ ea^/^ from (21) to bound (125) 
to be 



(126) 



There is also the field redefinition of vr in eq. (|27|) to take into account. In terms of C correlation 

)n 



functions this is equivalent to the field redefinition 

7f (") 



c 



{n} 



(127) 



We see that the second field redefinition is more important than the first by a factor a/e ^ 1, 
and so we concentrate on that. 

Let us see how the redefinition contributes to an n-point correlation function that we can 
denote by (C)red.: 

(nred = $^---$^(C^'"^^..C^'""^) . (128) 

mi rrin 

the fluctuations, we can conclude that this holds for all the orders, as we wanted to show. Explicit inclusion 
of spatial difF.s would modify the relations between C and tt to 



C^"^ = In 



a{t + T) 



a{t) 



/3 



(122) 



with . . . representing terms containing lower order terms in /3 and T . Eventiually we conclude that spatial 
diff.s can be safely neglected in the redefinition between Q and tt. 
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Since for the consistency of the effective theory we have a <^ C^^^"^, we see that the leading 
effect comes from substituting the lowest order possible field redefinition. Further, the struc- 
ture of the field redefinition with e ^ 1 and a ^ 1 suggests that the leading term is obtained 
by replacing only one of the C's with the quadratic field redefinition. We therefore have 

iOred ~ eosc«(C+') ■ (129) 

If n is odd, this is the leading effect; if n is even, we have to either go to higher order in the 
field redefinition or to compute the non-Gaussian correlation. Let us start with n odd. The 
signal to noise associated to this field redefinition is given by 

{C)red ^ /^2\l/2 ^ l/2/^2\l/2 

^^2^n/2 ~ ^osca(C ) ' ^ ea ' (C ) ' , (130) 

where we used that Cosc ^ ea~^/^. This is still less than taking the relation a <^ C^^^^ 
into account. In order to be detectable in the foreseeable future, we need this ratio to be 
larger than 10~^ (corresponding to the optimistic case of signal to noise capable of detecting 
Inl ~ 1)- By a coincidence this is of order (C^)^^^- Therefore the effect from the field 
redefinition is too small. 

If n is even, one way is to insert either a third order field redefinition, leading to a signal 
to noise of order 

1^ - eosca'iCY' < e«^/^(C^) < ^iCT' , (131) 

where in the second inequality we used the bound on a ^ C^^^^ coming from the consistency 
of the effective theory. This is again too small to be detectable. Alternatively we can use a 



non-Gaussian correlation function. From (37), this is schematically of the form 

(C)~eosc«'"''/'(C')"~S (132) 
which leads to a signal to noise of order 

1^ ~ eLa^-'^'iCr^' < ^'iCY' , (133) 

which is again too small to be detectable. We conclude that the effect of the field redefinition 
is too small to lead to a detectable signal. 

The effect of the multivertex diagrams and the collective effect of the lower order terms 
neglected above were discussed in [22]. However if the ratio of the resonance frequency with 
respect to the cutoff is as large as half, which is somewhat border line, then we can have 
up to the 10*'*-point function possible detectable, which is about the same value at which 
multi-vertex diagrams tend to become important 
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C Signal-to-Noise ratio for 3-point function 



To estimate the signal to noise ratio of the 3-point function we need to calculate 



S 



N 

yN 



V 



N 



d^ki d^k2 d^k3 (CkiCkaCka) (CkiCkaCks) 

(27r)3 (27r)3 (2vr)3 (Ck.Ck.CkaCk.Ck.Cka) 



(i^ki d^k2 d^ks 
(27r)3 (27r)3 {2n)^ 



(27r)V(ki + k2 + k3) 



AlBl{h){27rmO) 



-3 ' 



(134) 



where An and -B„ are given in (37) and where for simplicity we neglect the transfer functions 
and sky projection effects because we expect them not to affect the result significantly. Taking 
{2nf5^{0) = Vwe get 



((O) 
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2\n 



d^ki d^k2 d^ks ,3,3,, , , , Bl 

(2^(2^(2^(2^) ^^'^^+'^^+'^^)rri 



-3 • 



(135) 



Since the integrand depends only on the magnitude of /ci, ^2, ^3 we can use the delta function 
to eliminate angular variables of the last two integrals 



d^k2 d^ks 
(27r)3 (27r)3 



(27r)^5^(ki + k2 + k3 



3! 



k2dk2 



ki/2 



k2 



k^dk^ . 



(136) 



ki—k2 



Here we have also used the symmetry of the integrand with respect to interchange of ki,k2,k3, 
multiphed by 3!, and limited the integrals to the region ks < k2 < ki. 
Using X2,3 = A;2,3/^i, B^{ki) at leading order in a is given by 



1 



sin(aln(i^/fc^)) + 



^-2^-3 



(137) 



where 63 is a scale invariant function except for the sine term. Changing the variables of 

re obt 

d^ki 



integration from ^2,^3 to X2,a;3 in (135) we obtain 

2 



S_ 
N 



27r2(C2)3 



(27r)3 

d^kl 

(27r)3 



dX2 
1/2 

I{ki/K\ 



dxsX2X^b^{x2, xs, ki/k^ 



l-X-2 



:i38) 



where I{ki/k^) is the result of the integrals over X2 and X3 and in principle depends on ki. 
However, since the dependence on ki solely comes from sin^(Q;ln ((1 + X2 + x^) ki/k^)) in 63, 
for large a's it has already been approximately averaged to 1/2 by the X2 and 0:3 integrals. 
Therefore /(/ci/fc*) can be approximated to be a ki independent constant 



I{k^/K) ~ /(I) 



dxo 



1/2 



X2 



dX3 = 1/J 



(139) 



1-2:2 
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Similarly we can calculate the signal to noise ratio for the modifications of the 2-point 
function. Using 



B2ik) ~ — sin(aln(2A;/A;*)) 



(140) 



we have 



S 
N 



^2 V 



:V 



(27r)3 
(27r)3 



4sin2(a ln(2A;/A;,)) 



(141) 



where in the second line we have replaced sin^(a ln(2/c/fc=i<)) by its average. 

Therefore the signal to noise ratio in 3-point function compared to the modifications of 
the 2-point function becomes 



(#((C^))) _ v^/^/^(l)A3 v^a^e?'' 



(142) 



Even by taking a so high it reaches its upper bound (43) the ratio (142) reaches its maximum 
•s/3/2 ~ 0.87. For smaller a's, this number scales as a^/«saturation5 ^^"^ ^is long as the EFT is 
not strongly coupled it is expected to be much smaller than that. 

Similar techniques can be used to evaluate the signal to noise ratio for the 3-point functions 
of section 4.2.2 The results are plotted in Fig. |4l 



D Oscillating models and standard shapes 

Here we show that rapidly oscillating shapes (large a's) are almost perpendicular to the 
standard shapes and therefore the standard non-Gaussianty searches can not be used to 
constrain the model in question We trivially extend the definition of Cosine in reference 
[33] to incorporate the fact that oscillating shapes are not exactly scale invariant. The only 
difference is that we need to integrate over all three momenta in the observable range that 
here we take to be from kmm = 10~^ Mpc~^ to kmax = 10~^ Mpc~^. In [21] a different 
definition for Cosine was used which is better suited to the CMB data analysis, while the 
definition of [33] is more similar to the analysis from both the CMB and the Galaxy Survey 
points of view. In Fig. |5} we consider the shape 

1 k- 
sm(a\nK/k^,) -\ — cos(a\n K/k^)y , 
a kj 

3,1 ■'A 

^^We stress that the parametrically leadmg signal is in the 2-point function, and not in non-Gaussianities, 
that add very httle signal to noise. 
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COS(Fosc, ^local) COS(Fosc, ^'equ) 

0.04 r 



0.02 



-0.02 
-0.04 




Figure 5: Clockwise from top left we have plotted the cosine with local, equilateral and 
orthogonal shapes, in the valid range of a. We see that as a becomes large, the cosine goes 
to zero. 



and calculate the Cosine with Local, Orthogonal and Equilateral shapes as a function of a. 
Solid lines in the plot are showing the semi-analytic prediction for the cosine as a function of 
a. Indeed using the same analysis as in |20] it is easy to show that Cosine scales as 

TT 

Clocal ^ -^^^ ■ (143) 
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